Dynamical quantum phase transitions and broken-symmetry edges in the many-body 

eigenvalue spectrum 
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Many body models undergoing a quantum phase transition to a broken-symmetry phase that 
survives up to a critical temperature must possess, in the ordered phase, symmetric as well as 
non-symmetric eigenstates. We predict, and explicitly show in the fully-connected Ising model in 
a transverse field, that these two classes of eigenstates do not overlap in energy, and therefore 
that an energy edge exists separating low-energy symmetry-breaking eigenstates from high-energy 
symmetry-invariant ones. This energy is actually responsible, as we show, for the dynamical phase 
transition displayed by this model under a sudden large increase of the transverse field. A second 
situation we consider is the opposite, where the symmetry-breaking eigenstates are those in the 
high-energy sector of the spectrum, whereas the low-energy eigenstates are symmetric. In that 
case too a special energy must exist marking the boundary and leading to unexpected out-of- 
equilibrium dynamical behavior. An example is the fermonic repulsive Hubbard model Hamiltonian 
rl. Exploiting the trivial fact that the high energy spectrum of H is also the low energy one of —H, 
we conclude that the high energy eigenstates of the Hubbard model are superfluid. Simulating in a 
time-dependent Gutzwiller approximation the time evolution of a high energy BCS-like trial wave 
function, we show that a small superconducting order parameter will actually grow in spite of the 
repulsive nature of interaction. 

PACS numbers: 71.10.Fd, 05.30.Rt, 05.70.Ln 
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I. INTRODUCTION 

The temporal quantum evolution of an isolated macro- 
scopic system initially prepared in an out-of-equilibrium 
configuration is currently turning from an abstract con- 
cept, useful for discussing fundaments of quantum statis- 
tical mechanics, to a real phenomenon that can be ob- 
served and studied experimentally. This metamorphosis 
has been mainly driven by experiments on cold atomsj^ 
but it will be surely given further impulse in the near fu- 
ture by the fast progresses in time-resolved spectroscopy 
on condensed-matter systems. In fact, the early time 
dynamics (up to ~ 1 ns) of a material that is driven out- 
of-equilibrium e.g. by an intense ultra-short laser pulse is 
still uninfluenced by the environmental heat sink, hence 
it is to a good approximation the dynamics of an isolated 
system. 

An interesting class of experiments focuses on the dy- 
namics across phase transitions. In cold atom systems, 
this can be achieved by a sudden change of the exper- 
imental conditions, e.g. the depth of the optical lat- 
tice, which corresponds to suddenly altering the Hamil- 
tonian parameters, a unique opportunity offered by these 
systems. - In pump-probe experiments on real materials, 
one can instead tune the fluence of the pumping laser, i.e. 
the excess energy injected into the system.- That energy 
is supposed to undergo fast redistribution among all de- 
grees of freedom, first among the electrons (faster) , later 
among the phonons (slower), thus effectively heating the 
sample and raising its temperature. If the equilibrium 
phase diagram has a transition between a low tempera- 



ture phase and a high temperature one, the effective tem- 
perature rise could drive such a phase transition, though 
the system will eventually relax back to its initial equi- 
librium state by the coupling to the external thermostat. 
In reality, a dynamical phase transition in an isolated 
macroscopic system is not as trivial as one could imag- 
ine. Across a thermodynamic transition, ergodicity is ei- 
ther lost or recovered, hence it is not at all obvious that 
the unitary time-evolution of an initial non-equilibrium 
quantum state should bring about the same results as, at 
equilibrium, the adiabatic change of a coupling constant 
or of temperature. 

This is actually a central question at the basis of quan- 
tum statistical mechanics. Here, more modestly, we high- 
light a possible link between the ergodicity breakdown 
at a phase transition and the occurrence in the many- 
body eigenvalue spectrum of " broken-symmetry edges" , 
namely of special energies that mark the boundaries be- 
tween symmetry-breaking and symmetry-invariant eigen- 
states. 

Let us imagine a system described by a Hamiltonian %, 
which undergoes a quantum phase transition, the zero- 
temperature endpoint of a whole second-order critical 
line that separates a low-temperature broken-symmetry 
phase from a high-temperature symmetric one, see Fig.[T] 
where g is the coupling constant that drives the phase 
transition. For the sake of simplicity, and also because 
it will be relevant later, let us assume that the broken 
symmetry is a discrete Zi - all arguments below do not 
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FIG. 1. (Color online) Generic phase diagram of a model 
that possesses a symmetry broken phase in the temperature 
T vs. coupling constant g phase diagram, separated from a 
symmetric phase by a second order phase transition. 



depend on this specific choice - with order parameter 

i 

which is not a conserved quantity, i.e. a, TL ^ 0, and 

where V — >• oo is the volume and i labels lattice sites. 
Below the quantum critical point g < g c , Fig. [TJ the 
ground state is doubly degenerate and not Z% invariant. 
If | ty±) are the two ground states, they can be chosen 
such that 

(<F ± | a | *±) = ±TO, with to > 0. 

On the contrary, for g > g c , the ground state is unique 
and symmetric, i.e. the average of a vanishes. Since 
the symmetry breaking survives at finite temperature, 
see Fig. [TJ one must conclude that, besides the ground 
state, a whole macroscopic set of low energy states is 
Z2 not-invariant. The ergodicity breakdown in a sym- 
metry broken phase specifically implies that these states, 
in the example we are dealing with, are grouped into 
two subspaces that are mutually orthogonal in the ther- 
modynamic limit, one that can be chosen to include all 
eigenstates with (a) > 0, the other those with (a) < 0. 
Since the symmetry is recovered above a critical tem- 
perature, then there should exist a high energy subspace 
that includes symmetry invariant eigenstates. We ar- 
gue that there should be a special energy in the spec- 
trum, a "broken-symmetry edge" £*, such that all eigen- 
states with E < E* break the symmetry, while all eigen- 
states above E* are symmetric. In the case where the 
Hamiltonian has additional symmetries besides Z2, hence 
conserved quantities apart from energy, we claim that, 
within each subspace invariant under these further sym- 
metries, there must exist an edge above which symmetry 
is restored, even though its value may differ from one 
subspace to another as is the case in the model discussed 
in the next section. 

Let us for instance focus on any of these subspaces. 
The Z2 symmetry implies that all eigenstates are even or 



odd under Zi- The order parameter u is odd, hence its 
average value is strictly zero on any eigenstate, either odd 
or even. Nevertheless, we can formally define an order 
parameter e) of a given eigenstate | ^>e) through 
the positive square root of 



We denote by 
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Psb(E) 



(2) 



(3) 



the density of symmetry-breaking eigenstates, namely 
those with m(^> e) > 777,0, where too is a cut-off value 
that vanishes sufficiently fast as V — >• 00, being e = E/V 
and Ssb(c) their energy and entropy per unit volume. 
Seemingly, we define 



psiiE)=e VS SI (e) j 



(4) 



the density of the symmetry-invariant eigenstates, 
to($ ' e) < with Ssi(e) their entropy. We claim that 
there exists an energy E„ — Ve* that marks the mi- 
crocanonical continuous phase transition in that specific 
invariant subspace, such that 



lim Ssi(e) = 0, 
lim Ssb (e) = 0, 



for e < e» 
for e > e* 



(5) 
(6) 



We do not have a rigorous proof of the above statement, 
but just a plausible argument. Let us suppose to define 
the average m(E) > too of the order parameter over the 
symmetry-breaking eigenstates through 



777(e) = 



1 



Psb(E) 



rn(^ E >) S(E' - E) 



The actual microcanonical average is thus 



to(c) = 



Psb(E) 



Psb(E) + psi(E) 



777(e). 



(7) 



(8) 



In the thermodynamic limit V — > 00, hence mo — > 0, 
the continuous phase transition would imply the exis- 
tence of an energy e* such that, for e < e*, 777(e) ~ 
(e* — e)* 9 , where the exponent /3' may not coincide 
with the corresponding one in the canonical ensemble 
m(T) - (T c - Tf, while m(e > e*) = 0. Since the 
entropy ratio on the r.h.s. of Eq. ([8]) is either 1 or 
in the thermodynamic limit, we conclude that the criti- 
cal behavior comes from m(e < e») ~ (e* — e)^ , which, 
by continuity, implies m(e > e*) = 0, namely that there 
are no symmetry-breaking eigenstates with finite entropy 
density above e*, hence Eq. ©. This further suggests 
that symmetry-breaking and symmetry-invariant eigen- 
states exchange their role across the transition, which 
makes also Eq. (JS|) plausible. We do not exclude that 
symmetry-breaking eigenstates may survive above e*, or 
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vice versa for symmetric ones; we just state that, if they 
survive, their entropy is not extensive. 

We may also guess a generalization of the above pic- 
ture to the most common situation of a first order phase 
transition. In this case we expect two different edges, 
ei < e*. Below e\ the entropy density of symmetry- 
invariant states Ssi{c) vanishes in the thermodynamic 
limit, while above e*, the actual edge for symmetry 
restoration, it is Ssb(c) that goes to zero. 

If we accept the existence of such an energy thresh- 
old, then we are also able to justify, without invoking 
any thermalization hypothesis*^ why a material, whose 
equilibrium phase diagram is like that of Fig. [TJ may 
undergo a dynamical phase transition once supplied ini- 
tially with enough excess energy so as to push it above 

We mention once more that the above arguments are 
not at all a real proof. However, they can be explicitly 
proven in mean-field like models, like the fully connected 
Ising model that we discuss in section [TTJ There, we ex- 
plicitly demonstrate that the dynamical transition occurs 
because above a threshold energy there are simply no 
more broken-symmetry eigenstates in the spectrum. We 
believe this is important because it may happen that such 
an energy threshold, hence such a dynamical transition, 
exists also in models whose phase diagram is different 
from that of Fig. [TJ as we are going to discuss in section 

inn 



II. FIRST MODEL: THE FULLY CONNECTED 
ISING IN A TRANSVERSE FIELD 

We consider the Hamiltonian of an Ising model in a 
transverse field 

^E J q<^ q -^o\ 



(9) 



where N is the number of sites, 

< = E e " iqri < > 

i 

is the Fourier transform of the spin operators, and J q 
the Fourier transform of the exchange. In the (mean- 
field) fully-connected limit, J q = J<5 q o, the model ([9]) 
simplifies into 



U = -— (Tq^o - hero 



N 



S Z S Z -2hS x , 



(10) 



having set J = 1 and defined the total spin S = er /2. 
It turns out that the Hamiltonian Eq. (fTU| can be solved 
exactly. We shall closely follow the work by Bapst and 
Semerjian^i whose approach fits well our purposes. For 
reader's convenience we will repeat part of Bapst and 



Semerjian's calculations. We start by observing that the 
Hamiltonian (|10[) commutes with the total spin operator 
SS, with eigenvalue S(S+1), so that one can diagonalize 
T-L within each S S [0,N/2] sector, which contains 2S + 
1 distinct eigenvalues, each one g(S) times degenerate, 
where g(N/2) = 1 and, for S < N/2, 



9(S) 



N 
« + 5 



N 



N 
T 



1 



(11) 



which is the number of ways to couple N spin- 1/2 to 
obtain total spin S. We define 



s = -v(i-fc). 



(12) 



where k, for large N, becomes a continuous variable k 6 
[0,1/2]. For a given S, a generic eigenfunction can be 
written as 



(13) 



M=-S 



where | M) is eigenstate of S z with eigenvalue M £ 
[—5, S}. One readily find the eigenvalue equation^ 



E<f> E (M) = — — M 2 $ B (M) 



(14) 



y/S(S + 1) - M(M - 1) $ E (M - 1) 



+ y/S(S + 1) - M(M + 1) <&e{M + 1) 



We now assume JV large keeping k constant. We also 
define 



2M 



e [-l + 2fc,i-2fc], 



so that, at leading order in N, after setting E = Ne and 
$ E (M) = $ e (m), 

the Eq. (HU) reads 



e<I> e (m) 



-m 2 $ e (m) - - - 2k) 2 - m 2 



$ f m 

N 



$ E I m + — 



Following Ref0, we set 

$> e (m) oc exp [ — N cj> € (m)] , 



(15) 



(16) 



where the proportionality constant is the normalization, 
so that 



$ f I m ± — | oc exp 



N(j>Am± 



N 



<]> e (m) e 



=F2<(m) 
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Upon substituting the above expression into (fT5|) . the 
following equation follows 



f>' e (to) = - arg cosh 



hy/(l -2k) 2 



(17) 



For large N, in order for the wave function (fT6|) to be 
normalizable, we must impose that: (i) the 5Re0 £ (m) > 
0; (ii) the 3?e0 e (m) must have zeros, whicfo because of 
(i), are also minima. As showed in Ref. @, these two 
conditions imply that the allowed values of the energy 



iin(/_(m) 



< e < Max 



where 



/+(m) = -to 2 + h v/(l-2fc) 2 -TO 2 , 



/_(m) = -to 2 - h v/(l - 2k) 2 - m< 
At fixed fc, the lowest allowed energy is thus 



(18) 

(19) 
(20) 



0.5 - 



-0.5 




FIG. 2. The two function f+(m) and f-(m) for k = and 
ft = 0.9. The allowed values of the eigenvalues are those 
between the minimum of /_ and the maximum of /+. 



= min(/_(m)) = -(1 - 2fc) 2 - 



and occurs at 



h 2 

m 2 (fc) = (1 -2fc) 2 - — 



(21) 



(22) 



if h < h(k) = 2(1 — 2k), otherwise the minimum energy 
occurs at to = 0, 



1 (fc) = /_(0) = -/ l (l-2fc). 



(23) 



It follows that the actual ground state is always in the 
k = subspace and has energy 



-1 - x if h < h(0) = 2, 
-h if/i>/i(0). 



(24) 



In Fig. [2] we plot the two functions /+ (to) and /_ (to) 
for k = and /i = 0.9 < h(0). As shown by Bapst 
and Semerjian^ whenever /_(to) has a double mini- 
mum as in Fig. [2] any eigenstate with energy below 
e < e* = /_ (to = 0) is doubly degenerate in the thermo- 
dynamic limit N —> oo, being localized either at positive 
or at negative to, thus not invariant under Z^. On the 
contrary, the eigenvalues for e > e* are not degenerate 
and are ^2 symmetric. More specifically, any eigenfunc- 
tion $ £ (to) has evanescent tails that vanish exponentially 
with N in the regions where /_(m) > e and /+(to) < e. 
In Fig. [3] we show in the case e < e» the regions of 
evanescent waves. In this case, one can construct two 
eigenfunctions, each localized in a well, whose mutual 
overlap vanishes exponentially for N — > 00. This result 
also implies that the ground state, which lies in the k = 
subspace, spontaneously breaks Z 2 when h < h(0), hence 
h(0) — 2 — h c is the critical transverse field at which the 
quantum phase transition takes place. Such a degenerate 




FIG. 3. (Color online) For a given energy e < e„, we draw in 
red the regions e < /_ (m) where the wave function vanishes 
exponentially as N — > 00. 



ground state is actually a wave packet centered either at 
to = - h 2 /4 or at -■s/l - V/4, see Eq. ([22]) . 

More generally, it follows that for any given k and h < 
h(k) = 2(1 — 2k) there is indeed a "broken-symmetry 
edge" 



,(fc) = -/i(l-2fe), 



(25) 



that separates symmetry breaking eigenstates at e < e» 
from symmetric eigenstates at higher energies. In par- 
ticular, in the lowest energy subspace with k = 0, the 
edge is e*(0) = — h. Therefore, although in the simple 
mean- field like model Eq. (fTU|) . one can indeed prove the 
existence of energy edges that separate symmetry invari- 
ant from symmetry breaking eigenstates. We also note 
that subspaces corresponding to different k have different 
e*(fc), as we anticipated in the Introduction. 
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The role of the broken-symmetry edge in the 
quench dynamics 



N 

(S x ) = y sin. 



COS ( 



(S») 



N 



sm v sm < 




FIG. 4. (Color online) Pictorial view of the quench dynam- 
ics. The upper curve (red) corresponds to the Hamiltonian 
for t < characterized by hi = 0.8. We assume that for neg- 
ative time the wave function is the ground state for m > 0, 
which corresponds to the minimum of the curve at positive 
rrii, dashed vertical line. The intermediate curve (blue) corre- 
sponds to f-(m) for hf — 1.2, while the lowest curve (black) 
to hf = 1.5. The intercepts e = /_(mj) define the lower 
bound of allowed energies. 

The quench dynamics we examine corresponds to prop- 
agating the ground state at h — hi with a different trans- 
verse field h = hf > hi. In the specific case of a fully- 
connected model, this problem has been addressed by 
Refs. d and[ij In particular, it has been found^ that for 
hi < h c , i.e. starting from the broken-symmetry phase 
a dynamical transition occurs at hf = h* = (ft c + hi)/ 2 



For hf > ft*, the symmetry is dynamically restored 
while, below, it remains broken as in the initial state 
If, instead of a sudden increase from hi < h c to hf, om 
considers a linear ramp 



h(t) = 




fH)£, forte [0,t], 
for t > r, 



then the critical ft* increases and tends asymptotically 
to the equilibrium critical value ft c (0) for r — > oo,— as 
expected for an adiabatically slow switching rate. This 
result demonstrates that the dynamical transition is very 
much the same as the equilibrium one, and it occurs for 
lower fields h only because of anti-adiabatic effects, which 
is physically plausible. 

The dynamical transition can be easily discovered in 
the semiclassical limit. If we set 

N 

(s*) = - cose, 



for large TV it is safe to assume {S Z S Z } ~ (S Z )(S Z ) = 
N 2 cos 2 0/4, so that the total energy per spin, see Eq. 
(fill, reads 



N 



= — cos 9 — h sin 9 cos < 



(26) 



and is conserved in the unitary evolution. Through the 
Heisenberg equation iS z = [S z ,%] and upon expressing 
cos (j) as function of e and 9 by means of Eq. (|26|) , one 
finds the equation of motion of the order parameter m = 
cos 9 

d cos 9 



dt 



T2 y/(h sin 9 - e - cos 2 9) (ft sin 9 + e + cos 2 9) . 

(27) 

It follows that, until e < —h, the order parameter os- 
cillates around a finite value - the symmetry remains 
broken - while for e > — h it oscillates around zero - the 
symmetry is restored. The dynamical transition thus oc- 
curs when e* = —h, which we recognize to be the edge be- 
tween symmetry-breaking and symmetry-invariant eigen- 
states at k — previously defined. 

This correspondence can be shown to hold also in the 
exact solution of the previous section. To this end, imag- 
ine to start from the ground state at an initial hi < h c , 
which occurs in the subspace of k = 0, i.e. maximum 
total spin S = N/2, and let it evolve with the Hamil- 
tonian at a different hf > hi. We note that, since S 
is conserved, the time-evolved wave function will stay 
in the subspace k = 0. The ground state at hi is de- 
generate, and we choose the state with a positive av- 
erage of m, which is a wave-packet narrowly centered 
around m t = y/l — (hi/2) 2 . For very large N, when 
the contributions from the evanescent waves in the re- 
gions where SRe0 e (m) > can be safely neglected, the 
initial wave function decomposes in k = cigenstates 
of the final Hamiltonian with eigenvalues e such that 
f-(rrii) < e < /_)_(m,-). In Fig. 2] we show graphically 
the condition e > /_(mj) for hi = 0.8 and two values of 
hf = 1.2, 1.5. We note that the minimum value of the 
allowed energy for hf = 1.2 belongs to the subspace of 
symmetry broken states, while for hf = 1.5 it belongs to 
the subspace of symmetric states. It follows that, while 



for hf = 1.2 the long time average of m will stay fi- 
nite, for hf — 1.5 it will vanish instead. The critical 
hf = h„ is such that /_(to^) = /-(0) = —ft.*, namely 
ft* = 1 + fti/2 = (ft c + hi)/2, which is indeed the result 
of Ref. l9|. Therefore, the dynamical restoration of the 
symmetry is intimately connected to the existence of an 
energy threshold. When the initial wave function decom- 
poses into eigenstates of the final Hamiltonian that all 
have energies higher than that threshold, then the long 
time average of the order parameter vanishes although 
being initially finite. 
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We note that the dynamical transition in this particu- 
lar example is related to the equilibrium quantum phase 
transition, but it is actually unrelated to the transition 
at finite temperature^ In fact, at a given value of the 
transverse field h < h c , all eigenstates within the sub- 
spaces with k > h c — h are symmetric, while those with 
smaller k have still low-energy symmetry-breaking eigen- 
states. Since the degeneracy g (S 1 ), see Eq. (|TTj) . increases 
exponentially in N upon lowering S, hence raising fc, the 
cntropic contribution of the symmetric subspaces at large 
k will dominate the free energy and eventually drive the 
finite temperature phase transition. On the contrary, the 
quench-dynamics is constrained within the subspace at 
k = 0, hence it remains unaware that in other subspaces 
the eigenstates at the same energy are symmetric. 

This observation is important and makes one won- 
ders how the above result can survive beyond the fully- 
connected limit. Indeed, as soon as the Fourier transform 
of the exchange J q , see Eq. ©, acquires finite compo- 
nents at q ^ 0, states with different total spin, hence 
different k, start to be coupled one to each other - the 
total spin ceases to be a good quantum number, the only 
remaining one being the total momentum. Therefore, 
symmetry-breaking eigenstates at low k get coupled to 
symmetric eigenstates at large k. In this more general 
situation, there are to our knowledge no rigorous results 
apart from the pathological case of one-dimension, where 
the energy above which symmetry is restored actually 
coincides with the ground state energy, or, more rigor- 
ously, where excited states that breaks the symmetry do 
not have extensive entropy. However, we mention that 
a recent attempt to include small q ^ fluctuations on 
top of the results above, i.e. treating 

as a small perturbation of the bare Hamiltonian 

%o = --^ cfq<Jq - ha%, 

suggests that the dynamical transition in the quench does 
survive^ which we take as an indirect evidence that 
the energy edge does, too. We suspect the reason being 
that, once symmetry-breaking and symmetry-invariant 
eigenstates of same energy get coupled by (|28[). the new 
eigenstates, being linear combinations of the former ones, 
will all be symmetry-breaking. As a result, the broken- 
symmetry edge will end to coincide approximately with 
e*(fc = 0), i.e. with the maximum value among all the 
formerly independent subspaces. 

III. A DIFFERENT EXAMPLE: THE 
HUBBARD MODEL 

In the Introduction we inferred the existence of a 
threshold energy by the existence of a finite tempera- 
ture phase transition that ends at T = into a quantum 



critical point. However, we found in the previous section 
an energy edge above which symmetry is restored that is 
actually unrelated to the finite temperature phase tran- 
sition, while it is only linked to the quantum phase tran- 
sition. This suggests that such an edge could exists in 
a broader class of situations. In some simple cases, one 
can actually prove its existence without much effort. Let 
us consider for instance the fermionic Hubbard model in 
three dimensions, with Hamiltonian 

n = - E **i + H - c ) + U J2 "W> (29) 

where U > is the on-site repulsion. In the parameter 
space where magnetism can be discarded, the low energy 
part of the spectrum is that of a normal metal, hence 
all eigenstates are expected to be invariant under the 
symmetries of the Hamiltonian H, namely translations, 
spin-rotations and gauge transformations. We observe 
that the high-energy sector of the many-body spectrum 
obviously corresponds to the low energy spectrum of the 
Hamiltonian —%. The latter is characterized by an op- 
posite band dispersion but, more importantly, by an at- 
tractive rather than repulsive interaction. As a result of 
the Cooper instability, the low-energy spectrum of —% 
must comprise eigenstates | vp) with superconducting off- 
diagonal long range order, 

lim | d\ 4 d u d jf \ *) =| A | 2 > 0, (30) 

\t— j|— >oo 1 

which are not invariant under gauge transformations. 
Since these are identically the high-energy eigenstates of 
the original repulsive Hamiltonian Eq. (|29|) . it follows 
that the upper part of its many-body spectrum contains 
eigenstates with off-diagonal long range order. There 
must thus exist a special energy which separates low- 
energy gauge-symmetric eigenstates from high-energy su- 
perconducting ones. 

Suppose to prepare a wave function that initially has 
| A(i = 0) | 2 > 0, see Eq. §U$, and let it evolve with 
the Hamiltonian H, Eq. (|2"9"|) . If its energy is high 
enough, so that its overlap with the upper part of the 
spectrum is finite, then | A(t —> oo) | 2 > 0, otherwise 
| A(t — )• oo) | 2 — > 0, signaling once again a dynamical 
transition that, unlike the previous example of section HIl 
is now accompanied by the emergence rather than disap- 
pearance of long-range order. This surprising result was 
already conjectured by Rosch et al. in Ref. 11 as a possi- 
ble metastable state attained by initially preparing a high 
energy wave function with all sites either doubly occupied 
or empty at very large U. What we have shown here is 
that such a superfluid behavior is robust and it is merely 
a consequence of the high energy spectrum of H, which 
contains genuinely superconducting eigenstates. We fi- 
nally observe that these states have negative tempera- 
ture, hence are invisible in thermodynamics unless one 
could effectively invert the thermal population^ 

While the existence of a high-energy subspace of super- 
fluid eigenstates of the Hamiltonian (|2T)f is evident by the 



7 




500 1000 1500 2000 



FIG. 5. Time evolution of the superconducting order param- 
eter A(t) for a repulsive U = 0.2U C , where U c is the critical 
repulsion at the Mott transition, and four different initial val- 
ues Ao, the curves a, b, c and d. We also show in the inset 
the early time relaxation of A(t) in case a. Time is is units of 
the inverse of half the bandwidth. We observe that the curves 
b, c and d maintain a finite order parameter, unlike the curve 
a, although it corresponds to the largest initial Ao = 0.7. We 
also note that in the case d with the lowest Ao = O.f, the 
order parameter actually grows in time. 

above discussion, it is worth showing explicitly its conse- 
quences in the out-of-equilibrium dynamics. To this end, 
we prepare an initial wave function | ty(t = 0)) =| *&\), 
ground state of the BCS Hamiltonian 

Hbcs = + ]T Uj (c\ a c ja + H.c) +J2 (Ac^ + ff.c), 

i,j,<7 i 

(31) 

where the sign of the hopping is changed with respect to 
(f2U)) . and A is a control parameter that allows to span the 
alleged high-energy superconducting subspace of the re- 
pulsive Hubbard model by tuning the value of the initial 
order parameter 

A ee A(t = 0) = (* A | cj t <4 + c liCrt | * A ). (32) 

Since this model proved insoluble so far, we re- 
sort to an approximate method for simulating the 
unitary time-evolution: the time-dependent Gutzwiller 
approximation! 13 ' 14 In brief, the time-evolved wave func- 
tion is approximated by the form 

I *(*)> - II I *Aw)j (33) 

i 

where Vi(t) is a non-hermitian time-dependent varia- 
tional operator that acts on the Hilbert space at site i, 
and | v I / A(t)} the solution of the Schrcedinger equation 
with a BCS Hamiltonian like in Eq. pip but with time- 
dependent coupling constants tij and A. All variational 
parameters are determined through the saddle point of 
the action 

S = J '*(*(*) | -■«!*(*)>, (34) 
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I i I i I i I i I 
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\ 

FIG. 6. Long-time average of the order parameter, A, as 
function of its initial value, Ao. In the inset we show the 
energy, in units of half-bandwidth, with respect to the ground 
state one obtained within the Gutzwiller approximation, as 
function of Ao, as well as the broken-symmetry edge, the 
dashed line, corresponding to U = 0.2U C - We note that the 
energy, i.e. the average value of the Hamiltonian (|29[) on the 
ground state of (|3ip . is lower the greater Ao because of the 
sign change of the hopping in the two Hamiltonians. 



which is computed within the Gutzwiller approximation. 
The method has been described in detail elsewhere, see 
for instance Ref. Ila . and we just present the results here. 
For convenience, all calculations are done at half-filling, 
where they are much simpler, focusing on the paramag- 
netic sector, i.e. discarding magnetism, and assuming a 
flat density of states with half-bandwidth D that we take 
as the energy unit. However, in order to avoid spurious 
interference effects from the Mott localization that oc- 
curs above a critical U c , we concentrate on values of U 
safely below U c , in fact below the dynamical analogue of 
the Mott transition that is found when U > C/ c /23^ 

In fig. [5] we plot the time evolution of the order pa- 
rameter for U — Q.2U C and different initial values of A , 
see Eq. (|3"2"j) . as well as different energies, see inset of 
Fig. [BJ where we also show the time- averaged values of 
A(t) with respect to their initial values. We observe that 
at the lower energies, which actually correspond to the 
larger Ao, the order parameter, initially finite, relaxes 
rapidly to zero. On the contrary, above a threshold en- 
ergy, A(i) stays finite and even grows with respect to its 
initial value, see Fig. [6] If we identify that threshold en- 
ergy as our broken-symmetry edge, although rigorously 
it is only a lower bound, its dependence upon U is shown 
in Fig. H 

The fermionic Hubbard model that we have discussed 
so far it is only a very simple example where one can infer 
the existence of a threshold energy in the many-body 
spectrum that is not accompanied by any anomaly in the 
thermodynamics. In fact, the idea of comparing the low- 
energy spectrum of a Hamiltonian H with that of — 7i, 
which is actually the high-energy spectrum of Ti, is very 
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FIG. 7. Broken-symmetry edge measured with respect to the 
ground state energy for different values of U 



simple yet very effective. However, we must remark that 
model Hamiltonians %, like the Hubbard model above, 
are meant to describe low energy properties of complex 
physical systems. Therefore, it is not unlikely that the 
value of the threshold energies extracted by comparing 
the low energy spectra of rl and — rl could be above the 
limit of applicability of the model itself, in which case 
they would be devoid of physical relevance. 



possible because the eigenvalue spectrum of the sys- 
tem is characterized by a well defined energy, which 
we named broken-symmetry edge, that separates low 
energy symmetry-breaking eigenstates from high-energy 
symmetry-invariant ones. Once the initial energy exceeds 
such a threshold, a dynamical restoration of symme- 
try occurs that mimics the equilibrium phase transition 
as obtained by increasing temperature or adiabatically 
changing coupling constants. We explicitly demonstrated 
this mechanism in the fully-connected Ising model in a 
transverse field, a prototypical mean-field example of a 
quantum phase transition. 

After that, we realized that such edges might arise also 
when the equilibrium phase diagram does not display any 
order-disorder phase transition, and bring about unex- 
pected dynamical features. For instance, the repulsive 
Hubbard model has a high energy sector of genuinely 
supcrfluid eigenstates. By simulating in the Hubbard 
model the unitary evolution of a BCS wave function we 
surprisingly found that, if the energy exceeds a threshold, 
the initial superfluid order parameter not only survived 
in the long time limit, but may actually grew, despite the 
fact that the interaction is repulsive. 

In summary we believe that the concept of broken- 
symmetry edges that separate well distinct sectors of the 
eigenvalue spectrum is quite general and may be a useful 
framework of reference in the out-of-equilibrium quantum 
dynamics of isolated systems. 



IV. CONCLUSIONS 
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